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ABSTRACT 



We study pair production from a strong electric field in boost- 
invariant coordinates as a simple model for the central rapidity 
region of a heavy-ion collision. We derive and solve the renormal- 
ized equations for the time evolution of the mean electric field and 
current of the produced particles, when the field is taken to be a 
function only of the fluid proper time r = \Jt^ — z^. We find that 
a relativistic transport theory with a Schwinger source term mod- 
ified to take Pauli blocking (or Bose enhancement) into account 
gives a good description of the numerical solution to the field 
equations. We also compute the renormalized energy-momentum 
tensor of the produced particles and compare the effective pres- 
sure, energy and entropy density to that expected from hydrody- 
namic models of energy and momentum flow of the plasma. 



1 Introduction 



A popular theoretical picture of high-energy heavy-ion collisions begins with 
the creation of a flux tube containing a strong color electric field [Q]. The 
field energy is converted into particles as qq pairs and gluons are created 
by the Schwinger tunneling mechanism p|, ^, The transition from this 
quantum tunneling stage to a later hydrodynamic stage has previously been 
described phenomenologically using a kinetic theory model in which a rel- 
ativistic Boltzmann equation is coupled to a simple Schwinger source term 
[|, H, 1^, |]. Such a model requires justification, as does the use of Schwinger's 
formula in the case of an electric field which is changing rapidly because of 
screening by the produced particles. Our aim in this paper is to present a 
completely field-theoretic treatment of the electrodynamic initial- value prob- 
lem which exhibits the decay of the electric field and subsequent plasma 
oscillations. This approach allows direct calculation of the spectrum of pro- 
duced particles from first principles and comparison of the results with more 
phenomenological hydrodynamic models of the plasma. 

Although our approach is relevant to heavy-ion collisions at best only 
during the period when the produced partons can be treated as almost free, 
the details of hadronization are not expected to affect the average flow of en- 
ergy and momentum. Hence information obtained about energy flow in the 
weak coupling phase where our methods apply is translated to energy and 
momentum eventually deposited in the detector. Using further hadroniza- 
tion assumptions one can relate our results to the particle spectrum of the 
outgoing particles. 



Recently |T^, |T^, ^ ^ we have presented a practical renormalization 
scheme appropriate for initial-value and quantum back reaction problems in- 
volving the production of charged pairs of bosons or fermions by a strong 
electric field. In those papers the electric field is restricted to be spatially 
homogeneous, so that all physical quantities are functions of time alone. The 
method used to identify the divergences is to perform an adiabatic expan- 
sion of the equations of motion for the time-dependent mode functions. The 
divergence in the expectation value of the current comes from the first few 
terms in the adiabatic expansion and can be isolated and identified as the 
usual coupling-constant counter term. In this manner we were able to con- 
struct finite equations for the process of pair production from a spatially 
homogeneous electric field, and to consider the back reaction that this pair 
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production has on the time evolution of the electric field. 

In heavy-ion collisions one is clearly dealing with a situation that is not 
spatially homogeneous. However, particle production in the central rapidity 
region can be modeled as an inside-outside cascade which is symmetric under 
longitudinal boosts and thus produces a plateau in the particle rapidity dis- 
tributions 0, |lT], [1^, |1^]. This boost invariance also emerges dynamically 
in Landau's hydrodynamical model [14| and forms an essential kinematic in- 
gredient in the subsequent models of Cooper, Frye, and Schonberg [|1^ and of 
Bjorken [T^. The fiux-tube model of Low ||16| and Nussinov Jl^ incorporates 
this invariance naturally. 

Hence, kinematical considerations constrain the spatial inhomogeneity in 
the central rapidity region to a form that again allows an adiabatic expansion 
in a single variable, the fiuid proper time r = \Jt'^ — z^. In Landau's hydro- 
dynamical model one finds that after a short time Tq the energy-momentum 
tensor in a comoving frame is a function only of r. This is the essential 



assumption we make that allows us to apply the methods of ||18|, |19| ^ ^ 
to the heavy-ion collision problem. 

In our approach the initial conditions of the fields are specified at r = tq, 
that is, on a hyperbola of constant proper time. The comoving energy density 
is a function of r only.[| Then the electric field must also be a function solely of 
r. We can apply the adiabatic regularization method to identify and remove 
the divergences. The simplification introduced by the boost symmetry allows 
us to study the renormalization of this inhomogeneous initial value problem. 

A remarkable feature of Landau's model is the appearance of a scaling 
solution in which Vz = z/t. Alternatively, requiring invariance under longi- 
tudinal boosts |l^ also leads naturally to a scaling solution. In models 
based on Landau's ideas, one also assumes that it is possible to determine the 
particle spectra from the hydrodynamical fiow of energy-momentum by iden- 
tifying particle velocities with hydrodynamical velocities. Using our methods, 
we can assess the validity of these assumptions. We calculate numerically the 
evolution of the electric field in r and study the accompanying production 
of pairs. As it turns out, the Boltzmann equation, when modified to refiect 
quantum statistics correctly, does very well at reproducing the gross features 



"'^FoUowing the usage in hydrodynamics we shaU continue to refer to the coordinates t 
and 77 = ^ In (l^) fluid proper time and the fluid rapidity, respectively, or simply 

as comoving coordinates. 
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of the field-theoretic solution. This is not too surprising, since our field equa- 
tions are mean field equations (formally derived in the large- iV limit), and 
hence are semiclassical in nature. The field theoretic treatment presented 
here is the first term in a systematic expansion in the number of parton 
species or quark flavors. The next order in the expansion contains dynami- 
cal gauge fields as well as charged particles in the background classical field. 
Thus in the next order one can study equilibration due to scattering; one 
could also calculate, for example, lepton production and correlations in the 
evolving plasma from first principles. These systematic corrections require 
the use of the Schwinger-Keldysh formalism and will be discussed else- 
where. 

The outline of the paper is as follows. In Section 2 we formulate electrody- 
namics in the semiclassical limit in the (r, 77) coordinate system corresponding 
to the hydrodynamical scaling variables. This curvilinear coordinate system 
requires some formalism borrowed from the literature of quantum fields in 
curved spaces, which we review for the benefit of the reader unfamiliar with 
the subject. In Section 3 we perform the renormalization of the current 
using the adiabatic method of our previous papers. This is needed as the 
source term for a finite Maxwell equation of the particle back reaction on 
the electric field. Section 4 is devoted to the renormalization of the energy- 
momentum tensor of the produced pairs in the comoving frame. We also 
discuss there the relationship to the effective hydrodynamic point of view. 
The detailed results of numerical calculations in the (l-l-l)-dimensional case 
for both charged scalars and fermions are presented in Section 5. 

The paper contains three appendices. In Appendix A we present the 
necessary formulae used in computing the particle spectrum from the time 
evolution of the field modes. In Appendix B we prove that, for the boost- 
invariant kinematics of this problem, the distribution of particles in fluid 
rapidity is the same as the distribution of particles in particle rapidity. In 
Appendix C we reformulate the problem in the conformal time coordinate 
which turns out to be somewhat more convenient for actual numerical meth- 
ods. 



2 Electrodynamics in comoving coordinates 
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2.1 Scalars 



We consider first the electrodynamics of spin-0 bosons. We sliall use tfie 
metric convention ( — h ++) wliicli is commonly used in the curved-space 
literature. The action in general curvilinear coordinates with metric is 



(2.1) 



where 

V^(/) = {d^ - ieA^)(f) , F^, = d^A, - duA^ . (2.2) 

We use Greek indices for curvilinear coordinates and Latin indices for flat 
Minkowski coordinates. 

To express the boost invariance of the system it is useful to introduce the 
light-cone variables r and ?7, which will be identified later with fluid proper 
time and rapidity . These coordinates are defined in terms of the ordinary 
lab-frame Minkowski time t and coordinate along the beam direction z by 

z = Tsmhri , t = rcoshr7. (2.3) 

The Minkowski line element in these coordinates has the form 

ds^ = -dr^ + dx^ + dy"^ + r'^drf . (2.4) 

Hence the metric tensor is given by 

(7^, = diag(-l,l,l,r2). (2.5) 

with its inverse determined from g^'^g^p = 5p This metric is a special case of 
the Kasner metric [E3 . 



For our future use we introduce the vierbein V^f which transforms the 



curvilinear coordinates to Minkowski coordinates 



9,u = VXVai>. (2.6) 

where rjab = diag{ — 1, 1,1,1} is the flat Minkowski metric. A convenient 
choice of the vierbein for the metric ( |2.5| ) for our problem is 



i;» = diag{l,l,l,r} (2.7) 
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so that 

V;'^ = diag|l,l,l,^} . (2.8) 
Thus the determinant of the metric tensor is given by 

det V = = r . (2.9) 
The Klein-Gordon equation is 

^ V^{g'"'V^)Vu(P-m^(P = 0, (2.10) 



and Maxwell's equations read 

1 



where 



-MV^Fn=f, (2.11) 



j^' = C{-ie[(l)^d^^ - (9^0^)0] - 2eM^(0V)} • (2.12) 



Here C denotes the operation of charge symmetrization as discussed in p8| . 

We are interested in the case where the electric field is in the z direction 
and is a function of r only. In Minkowski coordinates the only nonvanishing 
components of the electromagnetic field tensor Fab are 

F,t = -Ft, = E{t) . (2.13) 

In the curvilinear coordinate system we have 

= (2.14) 

where we have chosen the gauge A,- = so that the only nonvanishing compo- 
nent of is Ari{T) = A. Using the relationship between the two coordinate 
systems we find that 

E(r) = ^ = --^. (2.15) 
r r dr 



In this gauge and coordinates the Klein-Gordon equation becomes 

-Or : 

T r 



+ -Jr - ^^(-9, - %eA{T)f - - dy' + ) 0(r, r],x,y) = 0. (2.16) 
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In order to remove first derivatives with respect to r, we define a rescaled 
field X by 

0=^X- (2.17) 
The Klein-Gordon equation for x is then 

[dl - ^ [(9, - zeA{T)f - ^] - - 9/ + m^) x{r, V,x,y) = 0. (2.18) 

We are interested in the solution of this field equation, where A is regarded as 
a classical field determined by the expectation value of the Maxwell equation 
(|2.11|) . This approximation ignores processes with photon propagators, and 
can be shown W^ to be the leading order in a large- iV expansion, where N 



is the number of flavors of the charged scalar field. 

These equations are to be solved subject to initial conditions at r = tq. 
We need to specify the initial value of the electric field and the density matrix 
describing the initial state of the charged scalar field. For the problem at 
hand it is sufficient to describe the charged scalar field by the particle-number 
density and pair-correlation density with respect to an adiabatic vacuum 
state [see ( |2.30|) below]. 

In the gauge we have chosen there is homogeneity in rj as well as in the 
directions x and y. This allows us to introduce a Fourier decomposition for 
the quantum field operator x proper time r. 



X(r, r^, xx) = / [dk] \MT)a^e^^-^ + fl^irp^e 



— ikx 



(2.19) 



where 



t^"^] - (2vr)3 ' 

k ■ X = kriT] + k_|_ ■ X_|_ , 

k_L = iK,ky), x± = {x,y). 



(2.20) 



The modes /k satisfy the equation 



dr^ 



(2.21) 
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with 



^lir) = 7rJ(r) + kl + m2 + ^, 

k„ — eA 
7r^(r) = . 



(2.22) 



We quantize the matter field by imposing commutation relations at equal r, 



d7' 



i6{ri — r]')^ (x_L — x'_i_) 



(2.23) 



Demanding that the usual commutation relations obtain for the creation and 
annihilation operators, 

[6k, bl] = al] = {2ny6\k - k') , (2.24) 

in d spatial dimensions, we find that /k must satisfy the condition 

/k(r)9./^(r) - /^(r)9./k(r) = t . (2.25) 

This latter condition is satisfied by a WKB-like parametrization of /k. 



/k 



(2.26) 



(2nk(r))V2 (2fik(r))i/2 ' 
Because of ( p.21|) , f2k must satisfy the same differential equation as appears 



m our previous papers 



18, 19, Eg 



Ifik 

2^]k 



3^k o 



(2.27) 



The dot denotes differentiation with respect to r. 

The only nontrivial Maxwell equation in = gauge is 



1 d 
r dr 



^ d .. . 
-— yl(r) 
r dr 



if) ■ 



(2.28) 



In terms of the charge densities A^^^ and and the correlation pair density 
F, we can write the Maxwell equation as 



dE 
1^ 



e / [rfk] ^ [1 + 2Ar(k) + 2F(k) cos(2yk)] • (2.29) 
J i Zk 
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The structure of ( |2.29| ) is similar to that of the equation found for the ho- 
mogeneous problem |19| . We have used the definitions 

(4,ak) = {27r)V{k-k')N+{k), 
(&^-k'&-k) = (27r)^5^(k-k')iV_(k), 
(&-k'ak) = (27r)'^5^(k - k')F(k) . (2.30) 

Note that we have taken A^+(k) = iV_(k) = A^(k) since the current compo- 
nent j"^ vanishes due to the Maxwell equation (Gauss's law), 

f = e[\^ [iV+(k) - iV_(k)] = -^d.iV^F^^) = . (2.31) 

J T \ —Q 



2.2 Fermions 

Let us now turn to the same problem in Dirac electrodynamics. The la- 
grangian density for QED in curvilinear coordinates (found, for example, in 
[E3]) gives rise to the action 



S = j d'^+^x (det V) Y^Y^^."^ + ^(Vj.^)7^^ - im^* - -F^^F'"' 

(2.32) 

where |]2^ 

\/^^^J = (^^ + T^-^eA^)^ (2.33) 
and the spin connection is given by 

r, = h^'vud.v: + t;X) , = ^b", 7I , (2.34) 

with V^^ the usual Christoffel symbol. We find that in our case (see P5|) 

= r, = ry = o 

r. = -^7V. (2.35) 

The coordinate dependent gamma matrices are obtained from the usual 
Dirac gamma matrices 7° via 

r = rv^'ix) . (2.36) 
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The coordinate independent Dirac matrices satisfy 



(2.37) 



From the action ( |2.32| ) we obtain the Heisenberg field equation for the 
fermions, 

(7^V^ + m) ^ = , (2.38) 



which takes the form 



^ = 0, 



(2.39) 



Variation of S with respect to yields the semiclassical Maxwell equation 



-9 



(2.40) 



If the electric field is in the z direction and a function of r only, we find that 
the only nontrivial Maxwell equation is 



1 dE{r) 
T dr 



2r 



^^7V* 



(2.41) 



We expand the fermion field in terms of Fourier modes at fixed proper 
time r, 

^(x) = J [dk] ^[6,(k)V'iL(^)e'''''e'P " + dl(-k)Ck.(^)e"''"'e-'P "]. (2.42) 



The ijj^g then obey 



d 



^is{r) = 0, 



(2.43) 



where vr.^ has been defined previously in (|2.22|) . The superscript ± refers to 
positive- or negative-energy solutions with respect to the adiabatic vacuum 



at r = tq. Following pOl, we square the Dirac equation by introducing 



d 



1 



Xs 



fks 



(2.44) 
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The spinors Xs are chosen to be eigenspinors of 7*^7^, 

7Vxs = KXs (2.45) 

with As = 1 for s = 1, 2 and = — 1 for s = 3, 4. They are normahzed, 

xlXs = 26rs . (2.46) 

The sets s = 1,2 and s = 3,4 are two different complete sets of hnearly 
independent solutions of the Dirac equation (see [^). Inserting ( p.44| ) into 
the Dirac equation ( p.43| ) we obtain the quadratic mode equation 

(^^ + ^^-^A,7r,j/±(r) = 0, (2.47) 

where now 

ujl = TTl + kl + m\ (2.48) 

If the canonical anti-commutation relations are imposed on the Fock space 
mode operators, then the ^p^g must obey the orthonormality relations 

<^t. - (2.49, 

T 

where r, s = 1,2 or 3,4. Using the orthonormality relations and ( p.44| ) we 
find, for a given k and s, 



00' 



/*"/^ + f*"!^ - iXsTT [f*''f'^ - /*"/^) = -5"^ (2.50) 



where = ± refer to the positive and negative energy solutions. Notice 
that from the off-diagonal relationship (a 7^ /?) we can express the negative- 
energy solutions in terms of the positive-energy ones. This fact we shall use 
repeatedly in what follows. 

Now we parametrize the positive-energy solutions f^^ in the same manner 



as in Eq. (3.1) of Ref. [|0 
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where f2ks obeys the real equation 



2n 



ks 



+ 



Ao Tin 



1 K 



2 1] 



ks 



A. 



r]^ 'ks 
ks 



^k(^) 



(2.52) 



Returning to the Maxwell equation and following steps (2.27)-(2.30) of 
Ref. 1201 we obtain 



1 dE{T) 
T dr 



^E/Mk](ki + m2)A.|/, 

s=l 



ksl ' 



(2.53) 



where we have taken the particle number A^(ks) and pair correlation density 
F(ks) defined by the analogs of eqs. ( |2.30| ) equal to zero for simplicity. 

Using the normalization conditions (|2.49|) -( ^.50|) we may express the 
mode functions and current in terms of the generalized frequency functions 

^ks) 



2|/, 



ksl 



ujI + fiL + 



^ks + As 71". 



2n 



ks 



2 AsTTy^r^ks 



(2.54) 



[See Eq. (3.7) of Ref. [20|. 



3 Renormalization 



3.1 Scalars 



The renormalization of the equations of the last section is straightforward, 
and is accomplished by analyzing the divergences in an adiabatic expansion 
of the differential equation for fik(T). We first present the regularization for 
the scalar case, where fik satisfies the differential equation ( p.27| ). The diver- 



gences of physical quantities such as the current and the energy-momentum 
tensor can be isolated by expanding Q in an adiabatic expansion. Up to 
second order this is given by 



1 1 f u 




(3.1) 
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The unrenormalized Maxwell equation is 



— r- 



dE 
dr 



e J [dk] ^ [1 + 2N{k) + 2F(k) cos(2?/k(r))] . (3.2) 



To study its renormalization in ci = 3 spatial dimensions we need to consider 
only the vacuum term, 



The adiabatic expansion (|3?T|) leads to 



T 



1 

— + 



+ 



(3.3) 



(3.4) 



The first term in ( p.4| ) is zero by reflection symmetry if we choose fixed inte- 
gration boundaries for the kinetic momentum fc^ — eA^j. The only divergent 
terms occur at second order. (Higher terms in the expansion have more 
powers of k in the denominator.) Using ( p. 221) and refiection symmetry, the 
right-hand side of (|3.4]) can be written as 



e / [dk] {k^ - eAY 



eA 



eA 



5eA 



\k, - eAf + 1 



4^6 



(3.5) 



Performing the fc^ and azimuthal angular integrations first we obtain 



487r2 



A2 



dk] 



A-^ 



+ 



A 



(3.6) 



where we have inserted a cutoff in the remaining transverse momentum in- 
tegration. The logarithmically divergent first term in ( p.6|) is 



247r2 
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In 1 + 



(3.7) 



We recognize the cutoff dependent infinite part as 

dE 



2 c- 2 

e be T 



dr 



(3.8) 
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where 

5e^ = (l/247r2) ln(A/m) (3.9) 

is the usual one-loop charge renormalization factor in scalar QED. Defining 
the renormalized charge via 

= g2^^ ^ e^5e^Y^ = e\l - eR^5e^) , (3.10) 

and using the Ward identity eE = enEn, we may absorb the divergence 
in the current into the left side of the Maxwell equation to obtain a finite 
renormalized equation suitable for numerical integration. 

In the d = 1 case there is no transverse momentum integration, and the 
charge renormalization is finite. The finite coefficient of the combination 
—A + A/t is r-dependent, and is given by 



127r ( + 



-1 



(3.11) 



The standard result in (1+1) dimensions is 6e^ = (127rm^)~^, and this is 
what is obtained for the spatially homogeneous problem by our method as 
well |jl9|. Thus, absorbing Se^ in the renormalization of the charge [see (|3.10|) ] 



leaves us with a (finite) r-dependent coefficient that is multiplied by the above 
combination, which appears now on both sides of the finite Maxwell equation, 
just as in the three dimensional case. This feature is peculiar to scalars 
in the r coordinate. The actual numerical solution of the scalar equations 
was performed in the conformal time coordinate u = In (mr) discussed in 
Appendix C. 



3.2 Fermions 



We turn to the renormalization problem in the spin-i case. The unrenormal- 
ized Maxwell equation is 



d_ 



m2)A, 



(3.12) 



Replacing VL and VL with uj and uj on the left-hand-side of ( |2.52D , we obtain 
the adiabatic expansion up to second order. 



TXTX 



ttMi 



vr 



+T — ^+--- (3.13) 



4 



4a;2 



2uj 
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Using this expansion in ( p.54| ) allows us to express the integrand of ( p.l2| ) in 
the form 



E(A:i + m2)(-2A.)) 



I/, 



ksl 



r 



(3.14) 

where Rk{T) falls faster than u ^ and so leads to a finite contribution to the 
current. Substituting ( p.l4| ) and using the definition ( |2.22| ) of tt^ and its first 
and second derivative into ( |3.12|) yields 



■A\dE 



A -2^ 

T , 



-e /[t/k]i?k(r) 



(3.15) 



where A is again the cutoff in the transverse momentum integral which has 
been reserved for last. 

The cutoff dependent term on the right-hand side is precisely the loga- 
rithmically divergent charge renormalization in 3 + 1 dimensions. Defining 
be^ = (l/Gvr^) ln(A/m) as usual and shifting this term to the left-hand side 
we obtain 

e— (1 + e^6e^) = -e^ / [dklRJT) , (3.16) 
dr J 

after multiplying both sides of the equation by e. The renormalized charge 
is 



1 + e^e^) 



(3.17) 



and the Ward identity assures us that crEr = eE. Hence we obtain the 
renormalized Maxwell equation 



dE 



R 



dr 



-en I [rfk]i?k(r) 



(3.18) 



where -Rk(T) is defined by Eq. ( p.l4[ ), and the integral is now completely 
convergent. 
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4 The Energy- Momentum Tensor and Effec- 
tive Hydrodynamics 

In our semiclassical calculations, we follow the evolution of the matter and 
electromagnetic fields. With these quantities in hand we can calculate, in 
addition to the particle spectrum, other physical quantities such as the en- 
ergy density and the longitudinal and transverse pressures. To obtain these 
quantities we derive the energy-momentum tensor in the comoving frame. 
We shall give explicit formulae here only for the fermionic case. 

The energy-momentum tensor for QED is obtained by varying the action 
in (PI). We find 



T: 



total 



-2 5S 



on _i_ rpem 



(4.1) 



with 



rpf ermion 
fiu 



rpem 

fJ,U 



« r 



pu 



(4.2) 



In the following we shall drop the superscript on the fermion part of the 
energy-momentum tensor where it causes no confusion to do so. 

We are interested in calculating the diagonal terms of the matter part of 
the energy-momentum tensor and in identifying them with the energy and 
pressure in the comoving frame. We begin with T^-, 



(0|T..|0) = ^E/Mk](A:l + 

' s=l •' 



m 



f*- — f- — f*+ A- f+ 

J ks i ks i ks J ks 



s=l ■ 



ar T 



(4.3) 



In the latter form only the positive frequency mode functions appear which 
is most useful for the adiabatic expansion below. Averaging over s = 1, 2 
and s = 3, 4 we may also write ( [4.3| ) in the form 



s=l •' 



(4.4) 



15 



This expression contains quartic and quadratic divergences in 3 + 1 dimen- 
sions present even in the complete absence of fields (the vacuum energy or 
cosmological constant terms) and a logarithmic divergence which is related 
to the charge renormalization of the last section. To isolate these divergent 
terms we express the integrand of T^-t as the sum of its second order adiabatic 
expansion and a remainder term, 

- 2E(^1 + m')^\Mf = + {kl + + Rrrik), (4.5) 

where RrrO^) falls off faster than u!~^ so that the integral over i?^T-(k) is 
finite. 

The first term in ( [4.5|) gives rise to a quartic divergence in 3 space 
dimensions (or a quadratic divergence in 1 space dimension) independent 
of the electric field and must be subtracted. The vr^ term in ([4.5|) gives 
rise to a quadratic divergence in 3 dimensions which must be likewise sub- 
tracted. However, in 1 space dimension this term yields a finite contribution 
to (0|TrT-|0) which must be retained, since it is r dependent, and cannot 
be absorbed into a cosmological constant counterterm. Strictly speaking, 
subtracting this term in 3 dimensions can only be justified by using a coordi- 
nate invariant regularization scheme for formally divergent quantities, such 
as dimensional regularization, where quartic and quadratic divergences are 
automatically excised. Only in such a scheme can one be certain that the 
divergence can be absorbed into a counterterm of the generally coordinate in- 
variant lagrangian in ( |2.32| ). The net result is that this term must be handled 
somewhat differently in the 3 and 1 dimensional cases. 

The term which is linear in A vanishes when integrated symmetrically. 
The term in (O) proportional to e'^A^ is logarithmically divergent in 3 di- 



mensions, and finite in 1 dimension. In fact it is precisely 

^^'^ ' (4-6) 

where is given by (|3.9| ) in 3 dimensions and by (127rm^)~^ in 1 dimension. 
In either case, it can be absorbed into a renormalization of the electric energy 
term of the stress tensor, 

A' 

2^2 



T;™ = — . (4.7) 
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by charge renormalization. When added to the electromagnetic term it gives 



1 + e^5e^ 



2r2 



El 



(4.8) 



Thus the exphcitly finite, renormahzed {T^-t) for the combined matter and 
electromagnetic system is 



(4.9) 



in three dimensions where /?,-,-( k) is defined in ([4.5| ). In one dimension the 
finite TT^ term that must be retained in the adiabatic expansion of ( |4.5| ) gives 
rise to an additional {127it'^)^^ on the right side of ([4.9|) . 
Turning now to T^^, the matter part is given by 



s=l 



(0|T,,|0) = 2r^ [dk]{kl + m^)Xs7r,\f^ 



(4.10) 



The adiabatic expansion of the integrand gives in this case 



+ |2 
ksl 



s=l 



+ 7r^r^i?k(T) . 



(4.11) 



Inserting this into ( [4.10| ) we find again that the quadratic and quartic di- 
vergences in 3 dimensions are independent of the electric field and that the 
terms proportional to A vanish, whereas the term proportional to A"^ again 
renormalizes the electromagnetic contribution to the energy-momentum ten- 
sor, 

rpem 

m - 2 ■ 

Therefore the fully renormalized total in 3 dimensions is 



(4.12) 



(T,,) = -]^EIt^ + t'J [dk] 7c,R^{t) . 



(4.13) 



As for (Ttt) there is a finite additional term in 1 dimension that must be 
added to this expression, which in this case is equal to a constant, (127r)~^. 
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For the transverse components of the matter stress tensor in 3 dimensions 
we have 



2r ,=1 



(4.14) 



In a precisely analogous manner one develops the adiabatic expansion of the 
integrand, isolates the divergences and the logarithmically divergent charge 
renormalization which combines with the electromagnetic stress, 



2 



A 

rpem _ rpem _ / h i cr\ 

^xx ^yy ' V^-^'^) 

and arrives at the renormalized form for the total (T^.^.) = (T^y) can be 
written as 

{Txx) = \\ W\j^-^\R,, - 7r,i?k] + \eI (4.16) 

where Rrr and i?k have been defined previously. This result for the transverse 
components may be obtained by consideration of the trace of the energy 
momentum tensor T^. From Eqs. ( [4.9|) , ([4.13|) , and ([4.16|) we note that this 
trace vanishes as m — ^> 0. 

Both the unrenormalized and renormalized total energy-momentum ten- 
sors are covariantly conserved, so that we have 

J- ~ I-' -^^Jmatter T ;i_i)em — u ; 

In the boost invariant proper time coordinates one may verify that this equa- 
tion takes the form, 



d^T^^+'^ + '^ = F,^f. (4.18) 



If we follow the standard practice and define the energy density and trans- 
verse and longitudinal pressures via, 

T^^ = diag(e,p_L,p±,P||r2), (4.19) 
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then the energy conservation equation takes the form 

d{eT 



dr 



+ P\\=Ej, (4.20) 



In most hydrodynamical models of particle production, one usually as- 
sumes that thermal equilibrium sets in and that there is an equation of state 
P|l = P||(e). For boost-invariant kinematics v = z/t, all the collective vari- 
ables are functions only of r and therefore p\\ is implicitly a function of e. 
In our field-theory model in the semiclassical limit there is no real scattering 
of partons and thus one does not have charged particles in equilibrium or a 
true equation of state. Hence the transverse and longitudinal pressures are 
different. In the next order in there is parton-parton scattering and it 
remains to be seen whether thermal equilibrium and isotropy of the pres- 
sures will emerge dynamically. Nevertheless, even in this order in one 
can define an effective hydrodynamics using the expectation value of the field 
theory's stress tensor ([4.19| ). Formally introducing the auxiliary quantities 



"temperature" and "entropy density" in a suggestive way in analogy with 
thermodynamics via, 

e + p\\=Ts] de = Tds (4.21) 
we find that the entropy density obeys, 

d{sT) _ Ejrj 



dr T 



(4.22) 



Notice that when the electric field goes to zero st becomes constant. If p\\ 
also goes to zero with r faster than 1/t , then er is constant when the electric 
field goes to zero. 

In standard phenomenological models of particle production such as Lan- 
dau's hydrodynamical model, one usually assumes that the hydrodynamics 
describes an isotropic perfect fluid whose energy momentum tensor in a co- 



moving frame has the form ( 4.19 ) with py = p±. Then one component of the 

(4.23) 



energy conservation equation becomes 

d{sT) 



dr 



which is of the same form as (|4.22|) in the absence of electric field, with the 



difference that the isotropic pressure enters into the definitions of the entropy 
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density and temperature of the fluid in ( [4.21| ). From these definitions one 
can also calculate directly the entropy density in the comoving frame as a 
function of r by 

Since we follow the microscopic degrees of freedom, we can also construct 
the Boltzmann entropy function in terms of the single particle distribution 
function in comoving phase space. This is done in Appendix A. 

Let us compare the energy spectra of the isotropic hydrodynamics with 
the results of our field theoretic approach. One quantity we wish to deter- 
mine is the amount of lab frame energy in a bin of fluid rapidity, dEiab/dr). 
For the isotropic hydrodynamics, in the lab frame one can write the energy 
momentum tensor in the form, 

T"* = (e + p)MV + pr7"'' (4.25) 

where = cosh?] and = sinhr/. Then calculating dEiab/dt] on a surface 
of constant proper time r we obtain, 

dEiab _ f rj^tadc^a 



drj J drj 

daa = A±{dz, —dt) = A_Lr/(cosh?7, — sinh?]) 



= A^e{Tf)rf cosh 7] , (4.26) 

where A± is a transverse size which in a flux tube model would be the trans- 
verse area of the chromoelectric flux tube. 

We show now that our microscopic hydrodynamics gives an identical re- 
sult, without any assumptions about thermal equilibrium. In fact, trans- 
forming the result of our fleld theory calculation ( |4.19|) to the lab frame, 

/ e cosh^ 7] + p|| sinh^ i] (e + py ) cosh i] sinh r] \ 
P± 

p± 

V (e + p\\ ) cosh r] sinh i] e sinh^ i] + cosh^ r] ) 

(4.27) 

and recalculating dEiab/drj again gives (|4.26|) , where e(r) is now explicitly 
calculable from the modes of the fleld theory. 



T 



•ab 
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In hydrodynamic models one assumes that hadronization does not effect 
the collective motion. If all the particles that are produced after hadroniza- 
tion are pions then the number of particles in a bin of rapidity should be just 
the energy in a bin of rapidity divided by the energy of a single pion namely, 



dN 1 dE, 



lab 



dr] rTT-TT cosh r] drj 

= a/-^ (4.28) 

To see if this formula is working in our parton domain we can instead use 
the mass of a parton in place of and check directly whether the spectrum 
of partons given by 

dN €(Tf)Tf 

— = ±112. 4.29 
dr] m 

agrees with explicit calculation of particle number in the field theory, as given 
in Appendix A. 

In order for the formula ( [4.29| ) to be independent of Tf we require that 
the electric field become vanishingly small and that the pressure go to zero 
faster than ^ at large r. Indeed we will find this is approximately true in 
the numerical simulations, the results of which we will present in the next 
section. 



5 Numerical results in (1 + 1) dimensions 

In this section we present the results of solving the back-reaction problem in 
two dimensions (proper time r and fluid rapidity t]) , and compare the results 
to a phenomenological Boltzmann-Vlasov model. In previous calculations 
using kinetic equations in flux tube models 0, |^, |^ it was assumed that the 
Schwinger source term (WKB formula) can be used by taking the electric 
field, hitherto constant, to be a function of proper time. However, in the 
Schwinger derivation the time parameter (which is not boost invariant) plays 
an implicit role. Therefore, it is not clear a priori if such a source term in 
the kinetic equations represents the correct rate of particle production. From 
the experience obtained in the spatially homogeneous case, we believe that 
if we know the correct source term, a phenomenological Boltzmann-Vlasov 
approach should agree with the semiclassical QED calculation. 
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The phenomenological Boltzmann-Vlasov equation in 3 + 1 dimensions 
can be written covariantly as 



Dt 



dpu 



Po9 



00 



dN 



(5.1) 



We shall write the transport equation in the comoving coordinates and their 
conjugate momenta, 



{'r,X,y,r]), = {Pr,Px,Py,Pr,) 



(5.2) 



In order to write the invariant source term in these coordinates, we begin 
with the WKB formula, which is 



dN 



[(-^)l/2cig0^3q]^2p^ 



±[l±2/(p,r)]e|E(r) 



X In 



1 ± exp 



7r(m2 + p5_) 
e\E(T)\ 



(5.3) 



if the (constant) electric field is in the z direction. The ± refers to the cases 
of charged bosons or fermions respectively. Our model for the spatially ho- 
mogeneous case consisted of applying this formula even for a time-dependent 
electric field. Here we will allow the electric field to be a function of the 
proper time, writing 

' (5.4) 



\E(t) 



In the spatially homogeneous case we assumed that particles are produced 
at rest, multiplying the WKB formula by S{pz). This longitudinal momentum 
dependence violates the Lorentz-boost symmetry. Here we assume ^ 
0] that the distribution is 5(p^), which is boost-invariant according to 
(p.7| ). Assuming boost-invariant initial conditions for /, invariance of the 



Boltzmann-Vlasov assures that the distribution function is a function only 
of the boost invariant variables (r, — y) or {r^p^). The kinetic equation 
reduces to 



±[l±2/(p,r)]er|E(r) 



X In 



1 ± exp 



e\EiT)\ 



(5.5) 
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Turning now to the Maxwell equation, we have from ( |2.28| ) that 

dE 



dr 



■ _ -cond _|_ -pol 
JV Jri "T Jri 



(5.6) 



where j™"*^ is the conduction current and is the polarization current 
due to pair creation |19|, The invariant phase-space in the comoving 
coordinates is 



1 d^p dp±dpri 



Thus in (1+1) dimensions we have 



(5.7) 



J 



cond 



■pol 
Jri 



2e 
2 



dpn 



dp^ ^Df 
p 



F^'l J 2'KTPr Dt 



, mer 



±[l±2f{p, = 0,r)] sign[E(r)]ln 

71 



1 ± exp 



vrm 



leEir) 



(5.8) 



Assuming that at r = Xj there are no particles, the solution of ( ^.5|) along 
the characteristic curves 



^ = eF (t) 
dr ' 



IS 



fip„T) = ± rdT'[l±2fip, = 0,T')]er'\EiT') 

•iTi 



(5.9) 



(5.10) 



X In 



1 ± exp 



7rm 



e\E{T') 

Thus the system ( ^.5D -( ^^ reduces to ^ 



5(p^-eA^(r') + eA^(r)). 



— r- 



dE 
dr 



dT'[l±2f{p, = 0,T')] 



Mr') - A{t) 



^[Air')^A{r)f + T- 



--T 



^ A related derivation (but without the generahty and fuhy covariance of the present 
one) for this model in terms of different variables can be found in j^, |^ . 
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X |^(r')|ln 



±[l±2/(p, = 0,r)] 




1 + exp 



TT 



Et) 



(5.11) 



In the above expression and in the following we introduce the dimensionless 
variables. 



eA-* A. 



eE_ 

TV? 



E. 



(5.12) 



For scalar particles we carry out the calculations in terms of the conformal 
proper time u (see Appendix C). This gives us more control over the physics 
at very early times, allowing us to choose a well-behaved initial adiabatic 
vacuum, corresponding to initial conditions 



(5.13) 



It is worth mentioning that the adiabatic vacuum in terms of u is not identical 
to the adiabatic vacuum in terms of r; they are related by a Bogolyubov 
transformation. The variable u is regular and improves numerical stability 
near the singular point r = 0. 

For an initial adiabatic vacuum state the renormalized Maxwell equation 
is given by 



dE 

du 



1 



^oo 27r 



1 



1 



(5.14) 



where 5e^ = (127rm^)^^. Equations ( |5.14| ) and ( |C7| ) define the numerical 
problem for the boson case. In solving ( ^.14| ) and (|C.71 ) we discretize the mo- 
mentum variable in a box with periodic boundary conditions, fc^ ±27m/L 
where L = 500 and n ranges from 1 to 3 x 10'*. The time step in u was taken 
to be 5 X 10"^ 

To compare the Boltzmann-Vlasov phenomenological model to the above 
semiclassical system, Eq. (|5.11|) is written in terms of the conformal proper 
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time variable u and becomes 



dE_ 

du 



du'[l±2f{p^ = 0,u')] 



A{u') - A[u) 



-2u' 



^[A{u') - A{u)f + 



X \E(u')\\n 



1 ± exp 



TT 



\Eiu' 



±[1 ± 2/(p, = 0, u)] .sign {E{u)) In 



1 ± exp 




(5.15) 



where dA/du 



In the fermion problem we perform the simulations in terms of r, and in 
this case the semiclassical problem is defined by Equation ( p.47| ) and by the 
Maxwell equation 



dE{T) 
dr 



2(4 



1 - el6e^ ^ 



2 r dk 



(5.16) 



In this 1 + 1 dimensions problem 6e'^ = (67rm^)^^, Ai = 1, and A2 = —1. The 
time step in r was taken to be 0.0005 with the momentum grid the same as 
for the scalar problem. 

Figs. 1 and 2 summarize the results of the numerical simulations of 
charged scalar particles in 1 + 1 dimensions. In Fig. 1 we show the time evo- 
lution of A{u), E{u) and ir^{u) for the case of E{u = —2) = 4 and /w? = 1. 
We see that in the first oscillation the electric field decays quite strongly, in 
contrast with the spatially homogeneous case 
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In the latter the degra- 
dation of the electric field results from particle production only. In this 
inhomogeneous problem our system expands, and hence the initial electro- 
magnetic energy density is reduced due to the particle production and this 
expansion. This degradation due to the expansion can be inferred by solving 
a classical system of particles and antiparticles that interact with a proper- 
time dependent electric field without a particle production source term. We 
also note that the larger the initial field, the smaller the period of oscillations. 

In Fig. 2 we compare the time evolution of the phenomenological model 
(dashed curve) with the result of the semiclassical calculation (solid curve). 
Initial conditions were fixed at u = — 2 [Figs. 2(a)-2(b)] and at m = 
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[Figs. 2(c)-2(d)]. We see that there is good agreement between the semi- 
classical solution and the Boltzmann-Vlasov model. This agreement also 
holds for different values of the initial electric field and for different coupling 
constants. The Bose enhancement increases the frequency, and hence a bet- 
ter agreement is achieved, as expected. We conclude that the WKB formula 
with Bose enhancement is a suitable source term for the boost-invariant prob- 
lem. It is worth mentioning that at very early times of the evolution (before 
r = 1) the particle production is negligible and the electric field falls very 
slowly, as can bee seen in Figs. 2(c)-2(d). 

In Figs. 3 through 9 we present the numerical results for fermions in 1 -|- 1 
dimensions. In Fig. 3 we compare the time evolution of the Boltzmann-Vlasov 
equation (dashed curve) with the results of the semiclassical calculation (solid 
curve), where the initial conditions were fixed at r = 1. All succeeding figures 
refer to these same initial conditions. In Fig. 4 we present the time evolution 
for re where e = T^-t-. We see that at large r this quantity oscillates around 
a fixed value. In Fig. 5 we show the time evolution of p/e where pr^ = T, 



In this lowest order calculation there is no true dissipation and no particular 
equation of state emerges from the time evolution, although there is some 
indication that p approaches zero faster than e. In Fig. 6 we present the 
evolution of the particle density dN/drj. After a short time (of order r = 15) 
the particle density reaches a plateau which doesn't change much in the 
subsequent evolution. This is consistent with the fact that the Schwinger 
particle creation mechanism turns off rapidly as the electric field decreases. In 
Fig. 7 we present the time evolution of re/ {dN/drj). One can see that at large 
r there is some indication that this ratio approaches the value of the mass 
(we choose m = 1), which agrees with the prediction of the hydrodynamic 
model discussed in Section 4. [see Eq. ( [4.29|) ]. This lends support to the 
idea that the pion spectrum can be calculated using ( |4.28| ). Defining the 
Boltzmann entropy by ( |A.10| ) of Appendix A, we plot ts as a function of 
r in Fig. 8. Notice that this quantity is roughly constant after r ~ 20, by 
which time particle production has nearly ceased. [Compare Fig. 6.] This 
agrees with the result expected from the hydrodynamic point of view, eq. 
( [4.23|) . Finally, in Fig. 9, we plot the effective "temperature," defined by the 
hydrodynamic relation ( |4.21| ), but using the Boltzmann entropy of Fig. 8. 

In conclusion, the present results using boost-invariant coordinates fall 
into line with previous studies [|19|, [2^] of boson and fermion pair production 
from an electric field with back-reaction. The renormalized field-theory cal- 
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culation is tractable, yielding oscillatory behavior for a relativistic plasma 
which can also be well described by means of a classical transport equation 
with a source term derived from the Schwinger mechanism modified by Bose 
enhancement or Pauli blocking. For the boost-invariant case the electric field 
decays much more rapidly than for cartesian coordinates, where the sole de- 
cay mechanism is transfer of energy to the produced pairs. The ability to 
use the transport-equation approximation for boost-invariant pair production 
justifies in part the use that has been made of this method of description in 
past studies of the production of the quark-gluon plasma, and opens the way 
for further applications in the future. 
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A The particle spectrum of fermions 



We present here the formulae for direct calculation of the fermion particle 



spectrum in the adiabatic method. For further discussion see |21[| . 

During particle production particle number is not conserved nor even 
uniquely defined. In terms of the adiabatic expansion of Sections 3 and 
4, however, a natural definition of an interpolating particle-number opera- 
tor suggests itself. One may simply expand the field in terms of the time- 
dependent creation and annihilation operators of the lowest order adiabatic 
vacuum, 



-ik-xi 



(A.l) 



where 



-1 



9i 



Xs^, (A.2) 



analogously to eq. ( 2.44 ) of the text, but in which the exact mode functions 
obeying ( |2.47|) are replaced by lowest order adiabatic mode functions 
g^g. The positive frequency adiabatic mode function is given explicitly by 
substituting uJk{T) for f2k(r) in the expressions (|2.51|) and (|2.54|) for f^^ in 
the text. 

The adiabatic basis functions and the corresponding Fock space particle 
annihilation and creation operators, as(k;r) and cj(— k;r) defined in this 
way are related to those defined in ( |2.42| ) by a time-dependent Bogoliubov 
transformation. This transformation is easily found by using the Dirac inner 
product, 

{u, v)= f dWu^^v = I rf'^xv^ . (A.3) 



Indeed by substituting the two expansions of the field operator in terms of 
the two different bases ( |2.42D and ( |A.1| ) we find 



a,(k; r) = (y^,(r)e^'^■^ ^) = a(ks; r)6,(k) + (3\ks- r)4(-k) 



(A.4) 



with 



a(ks;r) = 2{ 



dr dr 



9ks 



d_ 

dr 



ks 



s /ks} • 



(A.5) 
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Squaring this expression and using ( p.51| ) and its analog for the adiabatic 
mode function g^^, we arrive at 

|/?(k.;r)p = l-|a(k.;r)p 



1 - 

4|/il 



TT, 



2n 



ks 



(A.6) 

The expectation value of the number operator with respect to the adiabatic 
Fock space operators in ( [A .11 ) is then simply the sum over s = l,2ors = 3,4 
of 



N{ks,T) 



-- N+{ks)\a{ks; t)\^ + (1 - iV_(ks))|/5(ks; t)\^ 
+ 2Re{a(ks;r)/5(ks;r)F(ks)}. 
For initial conditions which correspond to the adiabatic vacuum, 

^ksiro) = u;k(i"o), 
^ks(To) = c^k(T-o) 

one can choose A^±(ks) = F{ks) = 0, thus N±{ks,T) = |/3(ks;T)p. 
the time-dependent particle number defined by ( |A.7| ) with ( |A.6|) has the 
property of starting at zero at r = if the initial state is the adiabatic 
vacuum state. At late times, when electric fields go to zero it approaches 
the usual out-state number operator. Thus we can identify the phase-space 
density as 

f{k„k^,r) ^ , , = mk„k^,T)\^ (A.9) 



(A.7) 



(A.8) 
Hence, 



dr] dkri dk±dx± 

The quantity defined in this way from first principles of the microscopic 
quantum theory agrees quite well (after coarse graining) with the single par- 
ticle distribution function /(p, r), obtained by solving the Boltzmann-Vlasov 
eq. From the single particle distribution one can construct the Boltz- 

mann entropy density in the comoving frame, 

^ ' " " /)ln(l-/)}. (A.IO) 



s{t) 



T 



dkn dk^ , ~ ~ 
l/ln/ 



{2n) 



1 



If / approaches a Fermi-Dirac equilibrium distribution at late r, then this 
Boltzmann entropy will agree with the quantity (|4.24| ) defined by the energy- 
momentum tensor in Section 4. 
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B Fluid rapidity distribution and particle ra- 
pidity distribution 

In hydrodynamical models one has a purely phenomenological description in 
terms of the collective variables — energy density, pressure and hydrodynamic 
four-velocity. Using a criterion for hadronization such as those described by 
Landau and by Cooper, Frye, and Schonberg one determines the 



particle spectra by making a further assumption that at break-up the fluid 
velocity is equal to the particle velocity. In our field theory treatment no 
such further assumption is needed as long as boost invariance holds and 
we determine the particle spectrum along a surface of constant r. In this 
appendix we will show the equivalence dN/drj = dN/dy. 

One ingredient in the proof is the fact that the transformation of coordi- 
nates to the (?7, r) system is a transformation to a local frame that moves with 
constant velocity tanh 7] with respect to the Minkowski center-of-mass frame, 
i.e., the comoving frame is not accelerated with respect to the Minkowski 
frame. Because of this, the total number of particles counted in that frame is 
the same as the number of particles in the Minkowski frame of reference. The 
second ingredient is that in order for the phase-space volume to be preserved 
under our coordinate transformation, we need to ensure that the transforma- 
tion in phase space is canonical in the classical sense of preserving Poisson 
brackets. 

Consider then the coordinate transformation 

/,2 2\l/2 -'-1 /t -\- Z 



2 \t-z. 

•Pr = EtjT — -pzlr prj = —Ez + tp. (B.l) 
The Poisson bracket is defined as 

SAB}-——-——-—— fB2l 
dp dx dE dt dp dx dt dE ' 

The Poisson brackets of these quantities are 

{r,r7} = 0, {Pr„Pr} = 0, 
K,r} = -1, {p„r/} = l. (B.3) 

We see that the above transformation is canonical. 
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The phase-space density of particles can be derived as shown in Appendix 
A, and it is found to be r^- independent. In order to obtain the rapidity 
distribution, we change variables from {rj, A;^) to {z, y), where y is the particle 
rapidity, 



Thus we have 



dN 



drjdk^ d\i±d'x± 
where the Jacobian is 



J 



dN 



dy dz (ik_|_(ix^ 



dk^/dy dkrf/dz 
drj /dy drj/ dz 



dkri drj 
dy dz 



(B.5) 



(B.6) 



Pr and can be rewritten as 

Pr = 

Pv = 



m±_ cosh(r7 — y) 
—Tm_[_ sinh(?7 — y) 



(B.7) 



The particle spectrum is calculated at a fixed value of r, soi] = sinh ^(z/r) 
rj{z). Thus functionally at fixed r we have 



Pr, + eAr, = kr, = kr,{r]{z) - y) . 
The chain rule then gives 



dkr. 



dz 



dkff drj dkrj drj 
drj dz dy dz 



(B.8) 



(B.9) 



At constant r, then, \J\ = dz/dkjj, which leads to the desired result 

dN _ dN 
dy drj 



(B.IO) 



Since the right-hand side of (|A.9| ) is rj independent, dN/drj is flat in rj. From 
( IB.IOI ) we conclude that the distribution dN/dy is fiat, as expected. 
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C Scalar electrodynamics in conformal coor- 
dinates 



Since r = is a singular point of our equations, we find it convenient to 
introduce the conformal time coordinate u viaQ 

mr = e" . (C.l) 

In (1+1) dimensions the line element reads 

2u 

ds" = -de + dz^ = —{-du" + dri") . (C.2) 

The transformation from the Minkowski t, z coordinates to the Kasner m, 77 
coordinates is conformal. We shall refer to u as the conformal proper time.0 
Instead of expanding the field x ^ ( |2.19[ ) we expand the field with 
the mode functions = fk/^/r, which satisfy 



^''17 I 1 1 1 1 1. ^ — /— /-I I ( II 

gk{r) = 0, (C.3) 



d'^Qk _^ I dgk _^ 



2 , (K - eAiT)y 
III I -|- 



r2 



dr"^ T dr 

where = k\ + m? . In terms of u the mode equation is 



d^gk 



du 



2 



+ wi{u)gk{u) = Q, (C.4) 



where 



2 

wl{u) = ^e^" + {k^ - eA{u)Y . (C.5) 
We parametrize gk in a WKB-like ansatz, 

gk{u) = ^ exp (-i f Wk{u')du'') , (C.6) 



•^In the radial Schrodingcr equation a singularity at r = prevents straightforward 
application of the WKB method. The transformation r — tq exp u maps the singularity 
from the origin to —00 and enables one to use the WKB approximation in terms of the 
new variable u p6| . Our situation is different, because in the vicinity of the singular point 
T = we can still use the adiabatic expansion in terms of r; the variable u is still helpful 
in avoiding numerical difficulties. 

For a general Kasner metric the conformal time is defined as r\coni = 
where g is the determinant of (|2.5|). 
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and again the real mode functions Wk satisfy the differential equation 



IWk _ 3 
2Wk AW, 



I + = wl{u) , (C.7) 



where the dot now denotes differentiation with respect to u. The Maxwell 
equation for an initial vacuum state is now 

or 

Performing the adiabatic expansion of Wk{u) for large k we find 



du^ du 



/ £A^ _ ^dAr\ ^ ^^.^^ ^ ^ 
\ du^ du J 

(C.IO) 



where 



in (3+1) dimensions, and 



5e2 = ^lnA/m (C.ll) 

= (C.12) 
127rm2 ^ ' 

in (1 + 1) dimensions, as expected. Note that the renormalization procedure 
introduces no additional interactions, in contrast to ( |3.11| ) in (77, r) coordi- 
nates. 
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FIGURE CAPTIONS 



FIGURE 1. Conformal proper time evolution of (a) the gauge field A{u), (b) 
electric field E{u), and (c) current jr]{u), for scalar particles in dimensionless 
units (|5.12| ). The initial conditions are that of adiabatic vacuum with respect 



to conformal u time ai u = —2 with initial electric field E{u = —2) =4.0 

FIGURE 2. Conformal proper time evolution of electric field (a) E{u), and 
(b) scalar particle current Jrf{u) with the same initial conditions as in Fig. 1 
(solid lines) compared to solution of the Boltzmann-Vlasov equation (dashed 
line), (c) and (d) are the same as (a) and (b) but for initial adiabatic vacuum 
conditions at m = 0. 



FIGURE 3. Proper time evolution of the system of (a) electric field E{t), 
and (b) fermion current Jr){T), for initial conditions at r = 1 with initial 
electric field E{t = 1) = 4.0 

FIGURE 4. Proper time evolution of Te(r) for fermions. 
FIGURE 5. Proper time evolution of p/e for fermions. 
FIGURE 6. Proper time evolution of dN/drj for fermions. 
FIGURE 7. Proper time evolution of rt/ {dN / drf) for fermions. 



FIGURE 8. Proper time evolution of Boltzmann entropy density s, multi- 
plied by r for fermions. 

FIGURE 9. Proper time evolution of the effective hydro dynamical "temper- 
ature" for fermions. 
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